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The role of turbulence in current generation and self-excitation of magnetic fields has been studied 
in the geometry of a mechanically-driven, spherical dynamo experiment, using a three dimensional 
numerical computation. A simple impeller model drives a flow which can generate a growing mag- 
netic field, depending upon the magnetic Reynolds number Rm = noaVa and the fluid Reynolds 
number Re = Va/v of the flow. For Re < 420 the flow is laminar and the dynamo transition is 
governed by a simple threshold in Rm > 100, above which a growing magnetic eigenmode is ob- 
served that is primarily of a dipole field tranverse to axis of symmetry of the fiow. In saturation 
the Lorentz force slows the flow such that the magnetic eigenmode becomes marginally stable. For 
Re > 420 and Rm ~ 100 the flow becomes turbulent and the dynamo eigenmode is suppressed. The 
mechanism of suppression is due to a combination of a time varying large-scale field and the presence 
of fluctuation driven currents (such as those predicted by the mean-field theory) which effectively 
enhance the magnetic diffusivity. For higher Rm a dynamo reappears, however the structure of the 
magnetic field is often different from the laminar dynamo; it is dominated by a dipolar magnetic fleld 
aligned with the axis of symmetry of the mean-flow which is apparently generated by fluctuation- 
driven currents. The magnitude and structure of the fluctuation-driven currents has been studied 
by applying a weak, axisymmetric seed magnetic field to laminar and turbulent flows. An Ohm's 
law analysis of the axisymmetric currents allows the fluctuation-driven currents to be identifled. 
The magnetic flelds generated by the fluctuations are significant: a dipole moment aligned with the 
symmetry axis of the mean-fiow is generated similar to those observed in the experiment, and both 
toroidal and poloidal fiux expulsion are observed. 

PACS numbers: 47.65.-l-a, 91.25.Cw 

Keywords: magnetohydrodynamics, MHD, dynamo, turbulence 



I. INTRODUCTION 

Astrophysical and geophysical magnetic fields are gen- 
erated by complex flows of plasmas or conducting flu- 
ids which convert gravitational potential, thermal, and 
rotational energies into magnetic energy [H, Q . A com- 
prehensive theory of the magnetohydrodynamic dynamo 
is elusive since the generating mechanism can vary dra- 
matically from one system to another. These variations 
arise from differences in free energy sources, conductiv- 
ity and viscosity of the conducting media, and geometry. 
Isolating and understanding the mechanisms by which 
self-generation occurs and the role of turbulence in the 
transition to a dynamo remain important problems. 

Dynamo action arises from the electromotive force 
(EMF) induced by the movement of an electrically- 
conducting medium through a magnetic field. This mo- 
tional EMF generates a magnetic field which, depending 
on the details of the motion, can either amplify or atten- 
uate the initial magnetic fleld. If the induced fleld rein- 
forces the initial magnetic fleld, then the positive feed- 
back leads to a growing magnetic fleld. The source of en- 



ergy for this dynamo is the kinetic energy of the moving 
fluid. The fluid may be driven by many different mecha- 
nisms such as thermal convection in a rotating body for 
the case of the Earth, or by impellers in liquid sodium 
dynamo experiments. 

When the system is turbulent, the turbulence likely 
plays an important role in the dynamo onset and the 
saturated state. In the saturated state, the backreac- 
tion of the self-generated magnetic fleld modifies the ve- 
locity field. It is well known that in hydrodynamics, 
turbulence converts large-scale motions into smaller and 
smaller eddies, a process known as a turbulent cascade. 
In magnetohydrodynamics (MHD), fluid turbulence can 
fold a large-scale magnetic field into smaller structures 
[3]. If the small-scale magnetic fluctuations are helical, 
they can, on average, generate a net EMF by interacting 
with the velocity fleld fluctuations and drive large-scale 
currents. When the magnetic fleld of this fluctuation 
driven current reinforces the original magnetic fleld self- 
excitation may be possible. Thus the generation of small- 
scale currents may explain observed large-scale magnetic 
flelds 
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Exact treatment of current generation in electrically- 
conducting fluids requires solving the MHD equations 
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governing the magnetic and velocity fields: 
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■V^B + VxvxB (1) 



J X B + pj/V^v -Vp + F(2) 



where p is the density, a is the conductivity, v is the vis- 
cosity, and p is the pressure. F is a driving term annotat- 
ing the sundary sources of free energy in the flow. Eqs.[T] 
and [2] are highly nonlinear and, without limiting assump- 
tions, are analytically intractable. Early dynamo theory 
focused on solving only Eq.[T]in the kinematic limit where 
the linear magnetic field stability of prescribed velocity 
fields was calculated to determine whether magnetic field 
growth was possible Due to advances in comput- 

ing power during the last decade, great progress has been 
made by performing numerical simulations of dynamos, 
which simultaneously solve the non-linear MHD equa- 
tions (Eqs. [T]and[21). These studies break into two sepa- 
rate classes: global simulations which attempt to model 
geophysical or astro phy sical dy namos such as the Earth 
and the Sun 0, [IH, M El Q [II , and simplified mod- 
els in which the geometry is simple enou gh t o uniquely 
identify particular physical effects [H ITtI Il8l| . 

The numerical simulations have been useful for study- 
ing magnetic field generation, even though they are still 
far away from being able to resolve the fluid turbulence 
of the actual systems. In particular, the role of the mag- 
netic Prandtl number, Pm — Rm/ Re, on threshold con- 
ditions for magnetic field growth is of importance for 
understanding magnetic field generation in the Earth, 
Sun, and in experiments. The linear self-excitation of 
the magnetic field is governed by the magnetic Reynolds 
number, Rm = fiQaLVo, where a is the molecular elec- 
trical conductivity, L is a characteristic size of the con- 
ducting region, and Vq is the peak speed. Hydrodynamic 
turbulence is governed by the fluid Reynolds, number 
Re = Vo£/i^, where £ is the characteristic size of the 
fiow. Simulations are capable of resolving the modest 
Rms needed to observe self-excitation, but not at the 
very high values of Re typical of low Pm dynamos. Re- 
cent studies in periodic boxes ^19,, 20] have focused on un- 
derstanding the generation of small-scale magnetic fields 
at low Pm, and simulations in cylindrical geometries with 
mean-flows have been performed [21] which show that the 
dynamo can be suppressed when turbulence is present. 
The periodic box simulations are particularly good at 
modeling inflnite, homogeneous turbulence, though these 
conditions are rarely, if ever, realized in actual astophys- 
ical or planetary contexts. Little work has been done to 
understand the dependence of large-scale magnetic fleld 
generation on Pm. 

To address more realistic models of astrophysical tur- 
bulence, research has turned to experiments. Experi- 
ments at Riga [H [H, [Hi and Karlsruhe [H, [11 [13 
use pumps to create flows of liquid metal through heli- 
cal pipes. These experiments are designed to be lami- 
nar kinematic dynamos, i.e. the average velocity fleld of 



the liquid metal is designed (through impeller and pipe 
geometry) to produce a magnetic fleld instability. The 
motivation for using liquid metal in the Riga and Karl- 
sruhe experiments is to allow helical flows, yet the con- 
duction and flow paths are not simply connected. Dy- 
namos in simply-connected geometries where the flow is 
unconstrained have yet to be demonstrated in an exper- 
iment. 

The self-excitation threshold of the Riga and Karlsruhe 
experiments is governed by the magnetic Reynolds num- 
ber. For particular flow geometries, the kinematic theory 
predicts a critical magnetic Reynolds number, Rmcrit, for 
self-excitation such that a dynamo transition is observed 
when Rm > Rmcrit- An important result from the Riga 
and Karlsruhe experiments is that the measured Rmcrit 
at which the dynamo action occurs is essentially gov- 
erned by the mean velocity field. Turbulence which was 
constrained by the characteristic size of the channel, £, 
apparently played little role. 

The kinematic theory does not provide a hydrodynam- 
ically consistent treatment of the fluid turbulence, and 
in simply-connected dynamo experiments the turbulent 
fluid motion will be pronounced. According to measure- 
ments in hydrodynamic experiments, the turbulent ve- 
locity fluctuations scale linearly with the mean velocity 
such that V = C {V). Mean field theory ^ predicts that 
turbulence can modify the effective conductivity of the 
liquid metal. Random advection creates a turbulent or 
anomalous resistivity governed by the spatial and tempo- 
ral scales of the random fiow. A reduction in conductivity 
due to turbulent fluctuations was observed at low mag- 
netic Reynolds number in liquid sodium [28| . The scaling 
of this turbulent resistivity is readily obtained by iterat- 
ing on the magnetic field in the nonlinearity of Eq. [1] 
and looking at the term that depends on gradients of B. 
For large Rm in a fluid with homogeneous, isotropic tur- 
bulence, the turbulent resistivity is proportional to v£^, 
and produces a turbulent modiflcation to the molecular 
conductivity. 
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where is a characteristic eddy size (presumed to be 
some fraction of L). The turbulent resistivity, as de- 
scribed above, operates even if there is no clear scale 
separation between the mean flow and the turbulence, or 
if mean quantities are nonzero. The turbulent conductiv- 
ity should be used for estimating the dynamo threshold: 
Rm = iiQaxVoL > Rmcrit results in a dynamo. Thus, 
the onset condition in a turbulent flow is governed by 



Rm > 
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1 - Ce^Rmcrit/L' 
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Note that the potentially singular denominator imposes 
a requirement on the effectiveness of a particular flow 
pattern for self-excitation; dynamos will only occur if 
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The small Pm of liquid metals implies large fluctuation 
levels and a turbulent conductivity. The influence of tur- 
bulent conductivity on self-excitation enters through the 
dimensionless number CRrricrit^v/ L. Through fluid con- 
straints, the flow-dependent parameters C, and Rrn^t 
can be manipulated. In the Karlsruhe experiment [29|, 
for example, £y is set by the pipe dimensions, rather than 
the device size hence iv/L can be taken to be a fraction 
of the ratio of the pipe dimensions to the device size. An 
upper bound would be iy/L = 0.06. We take C < 0.1, 
and Rrricrit ^ 40, hence CRrricrit^v/L < 0.24. We expect 
therefore that dynamo onset would be governed mainly 
by laminar predictions, as found experimentally. 

Turbulence plays a much greater role in governing self- 
excitation in geophysical and solar dynamos since there 
are no boundaries to keep small-scale flow from influ- 
encing the conducting region, and the values Pm in 
the Earth's core and in the convection zone of the Sun 
(Pm - 10-5 to 10-s and 10"^ respectively) 0, HH. 
This is also true for several experiments now underway 
which investigate magnetic field generation in more tur- 
bulent configurations [32, 33J. One such experiment, at 
the University of Wisconsin-Madison, uses two impellers 
in a 1 m diameter spherical vessel, to generate flows of 
liquid sodium with {V)> 15 ms"^. These^flows are pre- 
dicted by laminar theory to be dynamos [34j . The Madi- 
son experiment is expected to achieve Rm > 150 which 
exceeds Rnicrit by a factor of two. Such experiments 
have prompted a number of theoretical investigations 
into whether magnetic field generation is possible for the 
small Prandtl numbers of liquid metals in experiments 
without a mean flow [l^, [20] ■ The Madison Dynamo 
Experiment uses a simple two vortex flow which, accord- 
ing to a laminar kinematic theory, produces a transverse 
dipole magnetic field. The experiment presents a unique 
opportunity to test the numerical models; the spheri- 
cal geometry makes it particularly well-suited to being 
simulated, and the magnetic fields can be fully resolved, 
though the fluid turbulence, cannot be fully resolved by 
simulation since Re ^ 10^ in the experiment. 

In this paper, three-dimensional direct-numerical sim- 
ulations are used to model the dynamics of the exper- 
iment. The simulations are used to predict the behav- 
ior of the experiment and give guidance on what role 
turbulence might have on current generation and self- 
excitation. Section |TT] of the paper describes the numer- 
ical model used for solving the MHD equations. Sec- 
tion |TTT] describes results from Pm ^ 1 simulations where 
the flow is laminar and gives an overview of the large 
scale flow which is linearly unstable to magnetic eigen- 
mode growth. Section |IV] describes dynamos at lower 
Pm where the flow become turbulent. SectionlVlpresents 
simulations of a uniform magnetic field applied to axis of 
symmetry of the mean-flow in which turbulence gener- 
ated currents are investigated in subcritical flows. 



II. NUMERICAL MODEL 

The numerical model used in this paper solves the 
MHD equations in a spherical geometry, resolving the 
velocity fleld at the origin, and has a forcing term which 
mimics the impellers used in the experiment. The code 
is designed to simulate the behavior of a spherical liquid 
sodium experiment. Sodium, at 100 °C, is an electrically 
conducting fluid fully described by the imcompressible, 
resistive, viscous MHD equations. The code uses a spher- 
ical harmonic decomposition of the vector potential of the 
velocity and magnetic flelds in the 9 and cj) directions, and 
flnite difference representation in the radial direction. 

The dimensionless equations which govern fluid mo- 
mentum, magnetic induction, and solenoidal field con- 
straints are: 



— + i?TOo(v • V)v = -RmoVP (5) 
ot 

+ PmW'^v + RmoF + Rm^i x B, 

— = i?TOoV X V X B + V^B, (6) 

V-v = (7) 
V-B = 0. (8) 

In these equations, the time has been normalized to 
a characteristic resistive diffusion time of — ^qgo? 
where a is the radius of the sphere, and the velocity has 
been normalized to a characteristic velocity Vq so that 
i?mo = V^a^Qa. For the actual experimental device, the 
radius of sphere is a = 0.53 m; Rm^ — 100 corresponds 
to a characteristic speed of Vq « 15 ms~^. The vector 
field F is a stirring term of order 1 which models the im- 
pellers in the experiments. In practice, the velocity field 
resulting from the stirring term has a peak normalized ve- 
locity different from one. This resulting velocity field is 
used to define the resulting magnetic Reynolds numbers 
for a specific simulation, ie. Rm — max(V)Rmo. The 
relative importance of the magnetic and viscous dissipa- 
tion is expressed by Pm ~ v^-qCT which for liquid sodium 
is 10^5; the fluid Reynolds number Re is directly related 
to the magnetic Reynolds number by Re = Rm/Pm. In 
practice, the simulations have only been carried out for 
Pm > 0.1 which is sufficient to observed turbulence in 
the flows, but four orders of magnitude larger than in the 
experiments. 

Since the fluid is incompressible, the density evolution 
is unimportant and the pressure equation need not be 
evolved. Other numerical representations of a spherical 
MHD system solve for the pressure as a constraint on 
the flow [11] , especially in systems like stellar convection 
zones where compressibility is part of the dynamics [l0| . 
This simulation does not evolve the pressure explicitly; 
rather it solves for the vorticity. Taking the curl of Eq. 
[51 the expression for the time evolution of the vorticity 
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— RniaV X V X u: + RmoV x J x B (9) 

ot 



+ PmV'^uj + RmoV x F, and 
V = (Vx)"^^;. 



(10) 



The spectral decomposition is that of BuUard and Gell- 
man (BG), in which the velocity field is described by a 
spherical harmonic expansion of toroidal t and poloidal 
s fimctions Q, 

V = V X tf + V X V X sf, (11) 

and the magnetic field is described similarly 

B = V X Tf+ V X V X Sf, (12) 

where s, t, S and T are complex, scalar functions of 
r, 9 , and (j). This representation automatically satis- 
fies Eq. [71 To decompose Eqs. [11] [T^l each scalar func- 
tion is projected onto a spherical harmonic basis set, 
normalized by TV^,™ = ^(2£ + - ra)\I^A.-K{i + m)\: 
Yun{0,<l)) = NunPr{cos9)e'""t'. Ye,m is summed from 
m — 0, ...,£ and an extra factor of \/2 in Ni^m for m 7^ 
since the function represents a real field. The result for 
the magnetic field is 
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T{r,e,^,t) = Ni^rn}_^ ^ T,,„,(r,i)Pr(cos0)e 



im(p 



(13) 



00 7n—t 



S(r,0,0,i) = Ni,n}^ 2^ SLrn{r,t)Pp{co8e)e 



im<p 



(14) 



and similarly for the flow scalars, t and s. 

One advantage of the BG representation is that mul- 
tiple curls, which appear with every poloidal component 
of the vector flelds, reduce to Laplacians. The curl of 
a general solenoidal vector-field, W can also be repre- 
sented by two scalar functions of position, e and /. If 
W = Vxef+VxVx fr, then clearly 



V X W = V X {-V'^ff) + V X V X ef. 



(15) 



To determine the discretized version of the vorticity 
equations, Eq. [H] is expressed in terms of the toroidal- 
poloidal representation: 

a; = a;s+u;T = VxVxtf+Vx (-V^s)f. (16) 

By substituting this form into Eq. [9l the need to deter- 
mine boundary conditions on the vorticity is eliminated 
and only boundary conditions on the velocity field scalars 
are required. The evolution equations for the flow ad- 
vance become 

(9t 

— -PmWh = i?TOo [Gls + [V X F]s (17) 



dt 



PmVh = Rmo [G]^, + [V x F]^ , (18) 



where G signifles the sum of the advection and Lorentz 
forces. The fourth-order derivative can be computed by 
consecutive Laplacian operators. 

The Crank-Nicolson method is used to advance the lin- 
ear terms. This method implicitly averages the diffusive 
terms and computes a temporal derivative accurate to 
second order. The fluid advection term has a hyperbolic 
character due to the propagation of inertial waves making 
it advantageous to use an explicit advancement for non- 
linear terms. An explicit second-order Adams-Bashforth 
predictor-corrector scheme is used to advance the pseu- 
dospectral nonlinear terms. 

The pseudospectral method computes a function in 
real space and then decomposes it in spectral space. 
Pseudospectral methods avoid the complications of the 
full-spectral methods which rely on term-by-term inte- 
grations of spectral components (such as in the Galerkin 
method) and in general are much faster than full-spectral 
methods [36} . The pseudospectral method has the disad- 
vantage of introducing discretization error through alias- 
ing. This error is addressed by padding and truncating 
the spectrum [s^. 

The radial derivatives in the diffusive terms are com- 
puted through finite differencing on a nonuniform mesh. 
The finite difference coefficients for the and op- 
erators result in a nonsymmetric band diagonal matrix. 
The boundary conditions are folded into the matrix de- 
fined by the implicit linear operators with Gauss- Jordan 
reduction to ensure the matrix remains band-diagonal 
for ease of inversion. Using an optimized LU decompo- 
sition, the radial evolution is solved independently for 
each spectral harmonic. The scalar fields are then con- 
verted to real space and the nonlinear cross products are 
updated during predictor and corrector steps. 

The temporal evolution loops over a spectral harmonic 
index, thus individual boundary conditions for the re- 
spective harmonics are separately applied. The highest- 
order radial derivative in Eq. [TSj is fourth order, requir- 
ing four boundary conditions on the poloidal flow scalar. 
Since the velocity must permit a uniform flow through 
the origin, coordinate regularity implies 



ds(r = 0) 

s(r = 0), t(r = 0), — ^ = for ^ 7^ 1 



s(r = 0), t(r = 0), 



dr 

dhjr = 0) 
gj,2 



(19) 



= for £ ^ 1. 



For better numerical stability, the more stringent require- 
ment s, t — > is applied to turbulent simulations. The 
other boundary conditions are given by assumptions of a 
solid, no-slip boundary. For the poloidal flow, 



se,m{a) = 0, 



dr 



while for the toroidal flow 



0, 



t£,m(a) = 0. 



(20) 
(21) 

(22) 
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The discretization of the induction equation is straight- 
forward in hght of the method presented for the flow. Us- 
ing the magnetic field given by Eq.[T2l the induction term 
in Eq. [S] is projected into toroidal and poloidal compo- 
nents, grouping toroidal and poloidal contributions. The 
discretized expressions for the magnetic advance are: 



dT 

— - V^T = RmoNT 
ot 

— -V^S - Rm^Ns 
at 



(23) 
(24) 



where N is the spectral transform of the inductive term 
in the BG representation. Coordinate regularity gives 
the conditions for the magnetic scalar functions „ (r = 
0), T,,„(r = 0) = 0. 

The highest-order derivative of the magnetic advance 
is Oir^). Given the conditions on the magnetic field at 
the origin, a boundary condition on the magnetic field is 
needed at the wall. The outer surface of the Madison Dy- 
namo Experiment is stainless steel, modeled in the simu- 
lation as a solid insulating wall. The remaining boundary 
conditions are solved by matching the poloidal magnetic 
field to a vacuum field via a magnetostatic scalar poten- 
tial, and noting the toroidal field at the wall must be 
zero. This implies 



as,,™ , (^ + 1)^ _ 2£- 

1 ^^.-m — -T- 

ar a £ ^ 

T£,„i(a) = 0. 
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(25) 
(26) 



In Eq. [551 Ci^rn=0 if there are no currents in the sur- 
rounding medium, but can also be finite to represent a 
magnetic field applied by external sources. 

The timestepping, while unconditionally stable for 
the diffusive problem, is advectively-limited by an 
empirically-determined temporal resolution requirement 
of At < 5(Aa;)^ for a given spatial resolution. The 
spectral transform is the most computationally-intensive 
portion of the code requiring roughly 80% of the CPU 
time. The upper bound on the spatial resolution is: 
Ne ~ 64, ~ 128, Nr 400 which gives, with deahas- 
ing, £max = 42, or nearly 1000 modes. 

A forcing term, zero everywhere except at the location 
of the impellers in the experiment, drives the flow. The 
forcing term for the impeller model is 



F^{r,z) = p"^ sin {npb) + S 
Fz{r,z) = e sin (tt/jc) -f 7. 



(27) 
(28) 



The axial coordinate, z, and the cylindrical radius, p, are 
restricted to the region 0.25a < \z\ < 0.55a and p < 0.3 a. 
The impeller pitch, e, changes the ratio of toroidal (F^) 
to poloidal {Fz) force. The constants 5 and 7 control 
the axial force, and in this article are zero in all but the 
applied-field runs where stronger axial forcing is useful. 
The signs of F^, Fz are positive for z > and with F^ 
negative for z < creating the counter-rotation between 
the fiow cells. F is constant, which allows the input im- 
peller power F • V to vary. The region of the impellers and 



an example of the resulting flow are shown in Fig. [TJb) . 
These flows are topologically-similar to the ad hoc flows 
in several kinematic dynamo studies [1, [13, [111 , but differ 
in that they are hydrodynamically consistent solutions to 
the momentum equation. 



III. LAMINAR DYNAMOS 



The impellor model described above predicts dynamo 
action for sufficiently-strong forcing. For the particular 
case of Pm ~ 1, a laminar flow results, as shown in 
Fig. IHa). Starting from a stationary liquid metal, the 
evolution is observed to go through several phases. Ini- 
tially, the kinetic energy of the flow increases, as does 
the maximum speed (i?m) of the flow. The resulting 
Rm is above the critical value at which dynamo action 
is expected from kinematic theory. The magnetic field 
energy then increases exponentially with time. The mea- 
sured growth rate agrees with the growth rate A predicted 
by a kinematic eigenvalue code using the generated ve- 
locity fields and solving Eq. [5] for solutions of the form 
B oc iJoexp"^*. After this linear-growth phase, a back- 
reaction of the magnetic field on the flow is observed 
which leads to a saturation of the magnetic field. In this 
saturated state, the generated magnetic field is predomi- 
nantly a dipole oriented transverse to the symmetry axis, 
as seen in Fig. HJc). The m = 1 equatorially dominant 
structure of the dynamo (shown in Fig.[5]Jb)) is consistent 
with kinematic analysis. 

The orientation of the generated dipole is not con- 
strained by geometry and is observed to vary between 
simulations. When the saturation state is oscillating, (or 
damped with oscillations as shown in the Pm = 0.5 case 
in Fig. [5]) the dipole drifts around the equator and also 
undergoes 180° reversals. 

Self-excitation depends on the shape of the flow as well 
as the magnitude of Rm. An ideal ratio of poloidal to 
toroidal forcing exists (parameterized by e in Eq. I28p for 
which the critical magnetic Reynolds number is mini- 
mized as seen in Fig. [3] Minimizing Rmcrit makes the 
flow easier to attain experimentally. This optimal ratio 
can be understood from a simple frozen flux model de- 
scribing the stretch-twist-fold cycle of the dynamo (see 
Ref. [3^]). If the toroidal rotation is either too fast or 
too slow relative to the poloidal flow, the advected field 
is not folded back on to the initial field. 

For laminar flows, the backreaction is the result of two 
effects. First, an axisymmetric component of the Lorentz 
force is generated by the dynamo, slowing the flow and 
reducing Rm. Second, the flow geometry is changed such 
that the value of Rmcrit is increased. In saturation the 
growth rate is decreased to zero, as the confluence of Rm 
and Rmcrit in Fig. 2] shows. 
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FIG. 1: (a) A schematic of the Madison Dynamo Experiment. The sphere is 1 m in diameter and filled with 105-110 °C liquid 
sodium. High speed flows are created by two counter rotating impellers. Two sets of coils, one coaxial and one transverse to 
the drive shafts, are used to apply various magnetic fields for probing the experiment, (b) Contours of the toroidal velocity 
and contours of the poloidal flow stream function <I> where "Vpoi = V$ x V0, of the axisymmetric double vortex flow generated 
by the impeller model. The region of forcing is shown schematically along the symmetry axis. 



IV. TURBULENT DYNAMOS 

To investigate the effect of turbulence on the dy- 
namo transition, simulations are performed at lower Pm 
(higher Re). The flow changes from laminar to turbu- 
lent at Re « 420. Above this value, a hydrodynamic in- 
stability grows exponentially on approximately an eddy- 
turnover time scale Tc with a predominantly m = 1 spa- 
tial structure. Through nonlinear coupling, the instabil- 
ity quickly leads to strongly-fluctuating turbulent flows 
(a detailed discussion of the spectrum of the turbulence 
is deffered to Sec. V). Fluctuations about the mean flow 
exist at all scales, including variations in the large-scale 
flow responsible for the dynamo. The turbulence is inho- 
mogenous with boundary layers, localized forcing regions, 
and strong shear layers. 

The effect of these fluctuations on the dynamo onset 
conditions and on the resulting saturation mechanism de- 
pends sensitively upon the viscosity (parameterized by 
Pm) . Figure [5] shows an example of the broad range of 
dynamics exhibited by decreasing Pm, for an approxi- 
mately fixed value of Rm. The magnetic field dynamics 
fall into several regimes depending on Re: the laminar 
dynamo, a dynamo that starts turbulent, but relaminar- 
izes the saturated fiow, a turbulent dynamo, and finally 
a turbulent flow with no dynamo. At Pm ~ 1, the vis- 
cosity is large enough to keep both the magnetic field 
and velocity field fully laminar. The spectrum is dom- 
inated by the driven velocity field and by the magnetic 
eigenmode, and the saturation mechanism is the Lorentz 



braking and modification to the flow mentioned above. 

For Pm = 0.33, Fig.[5]shows a flow that is initially tur- 
bulent, but the saturated state is laminar. The turbulent 
saturation of the magnetic field results in a reduction in 
the fluctuations of the flow since the Lorentz braking has 
reduced Re below the hydrodynamic instability thresh- 
old (decreasing Re from 496 to 320). A hydrodynamic 
case, which evolves the fiow with B = 0, shows that flow 
turbulence persists without the addition of a magnetic 
field into the system. The Re threshold distinguishing 
the turbulent saturated state from a relaminarized satu- 
ration is Re ~ 630. 

If Rm is fixed near the experimental maxima while Re 
is increased beyond 700, no dynamo is observed. Despite 
the fact that the mean flow still satisfies the require- 
ments of a kinematic dynamo, the turbulent flow does 
not produce a growing magnetic field. Evidently, it is the 
turbulent fluctuations about the mean flow that prevent 
field growth. Using the mean flow (averaged over sev- 
eral resistive times) for the Pm = 0.22 (with Rm— 190, 
i?e=863) as a prescribed flow in a kinematic evolution 
of the induction equation gives Rmcrit ~ 93, as shown 
in Fig. [S] Even though the average flow has Rm well 
above Rmcrit there is no dynamo. However, when the 
conductivity is doubled such that Rm — 388 a turbulent 
dynamo reemerges in the simulation. Hence, an empiri- 
cal critical magnetic Reynolds number, Rmcrit,T can be 
defined which depends on Re through the degree of tur- 
bulent fluctuations in the flow. 

These results are consistent with the dynamo transi- 
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FIG. 2: (a) The kinetic and magnetic energy densities shown versus time with Rm =159 and Pm = 1. The time is in units of 
the resistive diffusion time = /loaa^. (b) The contributions to the total magnetic energy density from the m = 1 transverse 
dipole and the axisymmetric m — modes, (c) Magnetic field lines of a saturated dynamo state for a laminar flow with 
Rm = 150. 



tion being affected by the turbulent resistivity of Eq. 
[3l From analysis of the simulation results, the corre- 
lation time Tc, the eddy scale size £y, i and fluctuation 
levels C = v/Vq have been determined in order to es- 
timate the parameters in ar under the assumption that 
the homogeneous turbulence results roughly apply to this 
bounded, inhomogeneous flow. Typical volume averaged 
values measured in the Rm = 190, Re — 863 simula- 
tion are: = 0.022a, C = 0.45, and Tcorr- = 0.041tct, 
which yields a volume-averaged conductivity reduction 
of ctt/ct ~ 0.461. The diminished conductivity yields 
R'nicrit,T — 238. The results from all of the simulations 
are summarized in Fig. [7] which shows that an increasing 
i?m, at fixed Re reestablishes field growth where turbu- 
lent fluctuations had previously suppressed the dynamo. 
The dashed line in Fig. [7] shows that the correlation 
length and constant C increase with Re and eventually 
asymptote when the conductivity is effectively reduced 
by 70 %. 

The simulated turbulence has no de facto scale separa- 
tion. This might appear to pose a problem, given that our 
interpretation of the effect of turbulence is the introduc- 



tion of a turbulent resistivity, and the turbulent resistiv- 
ity of Mean Field Theory (MFT) is usually couched 
in scale separation arguments. However, it should be 
noted that the scale separation requirement associated 
with the a and (3 effects of MFT does not enter into 
the form of /3, but does guarantee that a < j3. This is 
because a is proportional to helicity whereas (3 is propor- 
tional to energy while a multiplies a lower derivative of 
the mean field than does (3. In this sense the lack of scale 
separation in the simulations is consistent with the ap- 
parent weakness of a turbulent-a effect in a regime with 
a turbulent resistivity. 

While the simulations are limited to Re < 2000 by 
computational speed and storage, we believe the sim- 
ulations capture the dominant effect since the fluctua- 
tions at the largest scales are the strongest contributers 
to the turbulent resistivity by the following argument. 
In Kolmogorov turbulence [40| the spectrum is E{k) oc 
e^/^fc~^/^, where e is the energy dissipation rate. Thus 
the turbulent resistivity goes as [J^^ q^'^E{q)dq]^^/'^ ^ 

e^/'^fep where ko is the wavenumber of the large scale 
eddies and fc^ is the dissipation scale wavenumber. In 
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FIG. 3: The dependance of the hnear growth rate of the least- 
damped magnetic eigenmode on impeUer pitch e. The transi- 
tion from damped to growing (A = point) defines the critical 
magnetic Reynolds number; Rm < Rnicrit to a growing mag- 
netic eigenmode. 
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FIG. 4: Rm and Rmcrit evolution during the saturation 
phase of a laminar dynamo. Rmcrit is calculated from lin- 
ear stability for each instantaneous velocity field during the 
simulation. In saturation, Rm = Rmcrit- 



K41 turbulence, fc^ oc Re^^'^, as Re becomes large in 
comparison to Rm, the effect of turbulent fluctuations on 
conductivity will asymptote to a fixed value. It should be 
noted that the simple dimensional analysis used for esti- 
mating the turbulent resistivity reflects isotropic homoge- 
nous turbulence and is derived in the limit that there is 
no mean flow; this dynamo relies almost entirely on the 
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FIG. 5: The magnetic and kinetic energy densities for runs 
with fixed Rm {Rm = 165 ± 3%) but different Pm versus 
time in Tc- Note that Pm = 0.33 shows a relaminarization 
of turbulent flow while Pm = 0.22 is barely amplifying the 
initial noise and is shown multiplied by 50,000. 



(b)3 


















/ Dynamo 


Subcritical/ 



2 3 4 5 
Time Ix^j] 



70 



80 90 100 110 
Rm 



FIG. 6: (a) Rm as a function of time. Taking a series of 
flows over the range shown the mean flow is calculated giving 
Rm = 193 and Re = 863. (b) The kinematic analysis of the 
average flow, where only Eq. [6]is evolved, and Rm is varied, 
to determine Rmcrit = 91.5. 
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FIG. 7: Re — Rm phase diagram. A number of simula- 
tions whose hydrodynamic and final saturated states are doc- 
umented in Fig. [5l Rmcrit for the mean fiow (V) is essen- 
tially independent of Re, while the effective dynamo thresh- 
old grows with Re. The dashed line shows the qualitative 
behavior of the dynamo threshold in turbulent fiows (Vi). 
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FIG. 8: The growth rate of the dominant eigenmode calcu- 
lated for a time-series of flow profiles. 



presence of a mean flow. 

An alternative viewpoint, consistent with the phe- 
nomenological interpretation of enhanced resistivity put 
forward here, is that the large-scale variations in the ve- 
locity field are continuously changing the spatial struc- 
ture and growth rates of the magnetic eigenmodes of the 
system. A more thorough treatment of the dynamic vari- 
ation of dominant eigenmodes can be found for a slightly 
different problem in [4l|. Two effects can be important. 
First, the instantaneous growth rate of the least damped 
eigenmode fluctuates between growing and damped. For 
a run with Rm — 193, Re — 893 and e — 0.4, shown 
in Fig. [51 a dynamo occurs only when the flow spends 
sufficient time in phases which are kinematic dynamos. 
The kinematic growth rate is most often positive, con- 
sistent with the time averaged flow having growing mag- 
netic field solutions, yet the modifications made to the 
flow during the subcritical periods are sufficient to stop 
the dynamo. Second, the turbulence couples energy from 
the growing magnetic eigenmode into spatially-similar 
damped eigenmodes. As the flows evolve, the spatial 
structure of the eigenmodes change. The magnetic field 
structure of a single eigenmode at some previous instant 
in time must be described in terms of several modes 
after the flow changes. This transfer of energy from 
the primary mode is equivalent to enhanced dissipation. 
Analysis of the eigenmode structure shows that the least 
damped eigenmode during a nondynamo phase in Fig. 
|5] varies between the marginally stable nonaxisymmetric 
dipole and a stable axisymmetric magnetic mode. Thus 
the flow imparts energy to a spatially similar, but distinct 
magnetic eigenmode. 

Finally, it should be noted that distinguishing between 
growing and damped magnetic fields is difficult in the tur- 
bulent simulations. Typically, the turbulent runs have 
been limited to durations of less than 10 To-. The tran- 
sition may also be considerably more complex as seen in 
Fig. [5] where magnetic energy of the Pm = 0.22 simula- 



tion may show intermittent growth near Rrricrit.T- The 
simulations are thus consistent with intermittent excita- 
tion of the dynamo eigenmode by the mean fiow. The 
peak magnetic energy is limited by the magnitude of the 
initialized noise the simulation is started with instead of 
the backreaction with the flow. This effect is especially 
relevant when the magnetic field is sustained by an ex- 
ternal source as shown in Section IVl 

A dynamo still occurs in these flows for sufficiently 
large Rm (keeping Re fixed, growing magnetic energy is 
detected for sufficiently large Pm) . An example of a time 
evolution and spatial structure of a saturated turbulent 
magnetic field is shown in Fig. [5] for a simulation with 
Rm — 337 and Re = 674. An m = 1 transverse dipole 
field is still present, as in the case of the laminar dynamo 
in Fig. [21 however this turbulent dynamo is now domi- 
nated by the presence of a large m = field, aligned with 
the axis of symmetry of the impellors as shown in Fig. 
Mjo). This component, by itself, would appear to violate 
Cowling's theorum [42|, and so it must result from non- 
axisymmetric components of the velocity field and the 
magnetic field. Thus it appears probable that the na- 
ture of dynamo has changed fundamentally from a simple 
eigenmode driven by the two vortex flow, to a dynamo in 
which the turbulent fluctuations may be responsible for 
generating the large-scale magnetic field. 



V. SIMULATIONS OF A SUBCRITICAL 
TURBULENT FLOW WITH A WEAK, 
EXTERNALLY-APPLIED MAGNETIC FIELD 

As a means of further highlighting the different physics 
and conditions underlying turbulent and laminar dy- 
namos, subcritical flows arc simulated with focus on the 
potential role of fluctuation driven currents in the self- 
excitation process. Subcritical flows have Rm < Rmcrit, 
and are not expected to lead to self-excited magnetic 
fields. The MHD behavior is investigated by applying 
a magnetic field which is generated by currents flowing 
in coils external to the sphere. The configuration studied 
is similar to the set of experiments described in Ref. [i^ , 
and is deliberately set up as an axisymmetric system in 
which fluctuation driven currents can be easily detected. 

The numerical technique employed is similar to the 
dynamo simulations described above in all but one re- 
spect, namely a different boundary condition is used 
with Ce^m 7^ in Eq. [23 These boundary conditions 
match the magnetic field to a scalar magnetic potential 
B — — V^m, which is valid in the region between the sur- 
face of the sphere and the external magnets. <& satisfies 
Laplace's equation and its solution is well known: 

$,„(r, e,c^)^J2 (^^.™^' + DL,nr-^'+^^) Y^e, 0), 

(29) 

where Y^{9, (j>ys are the spherical harmonics. The I?^,m 
terms represent the magnetic field generated by currents 
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FIG. 9: (a) The energy density of a turbulent dynamo with Rm — 337 and Re — 674. (b) The energy density of the 
axisymmetric magnetic field (m — 0) and the nonaxisymmetric dynamo (m = 1). (c) Magnetic field lines of the turbulent 
saturated dynamo 



in the sphere, and the coefficients Ai^m can be chosen 
to describe a magnetic field of arbitrary shape and ori- 
entation applied by currents external to the sphere. In 
this paper and in the simulations described below, a uni- 
form magnetic field is applied along the symmetry axis 
of the forcing terms, and is characterized by a single co- 
efficient Aifi, all higher order terms being zero. The 
appHed magnetic field, i?i,o, is weak enough so that it 
does not alter the large-scale flow. The strength of the 
applied magnetic field is moderated by keeping the Stu- 
art number N = aaB^ q/ pvQ < 0.1. In sodium, with 
a Rm ~ 100, N ~ 0.1 would correspond to an applied 
field of 156 gauss. The applied field for these simula- 
tions is uniform and applied along the impeller axis with 
Bifi ~ 57 gauss and N ~ 10^^. However, since the ve- 
locity fluctuations decrease with scale, the Stuart num- 
ber increases with scale, indicating the Lorentz force may 
influence small-scale fluid motion. Examples of such sim- 
ulations are shown in Fig. I10[ where the kinetic energy 
and magnetic energy are shown for laminar and turbulent 
runs. 

For the laminar flow, the induced currents and result- 
ing magnetic field are purely due to the magnetic field 
interacting with the mean flow, as seen in Fig. Illf a). 



Two main effects are observed. First, induced toroidal 
currents compress lines of poloidal magnetic field near 
the axis of the device. The lines are pulled outward at 
the poles and inward at the equator. The net result is 
a reduction of the poloidal field strength at the equator 
in the outer region, and a large amplification at the axis 
(the peak poloidal field is 18 times the applied field). Sec- 
ond, poloidal currents generate a toroidal magnetic field. 
These currents are generated by the well-known omega 
effect of dynamo theory whereby differential toroidal ro- 
tation of the fluid is able to stretch the field into the 
toroidal direction [l| . The amplitude of the peak toroidal 
field is greater than 6 times the applied field. 



The transition to turbulence is still characterized by 
the same Re ^ 420 threshold described above, since the 
Stuart number for the applied magnetic field is small. 
Below this threshold, the nonaxisymmetric part of the 
flows is negligible while above this threshold nonaxisym- 
metric fluctuations in both B and V can be as large as 
40% of the mean values. The geometry of the simula- 
tions (axisymmetric drive terms aligned with the applied 
magnetic field) makes it possible to separate mean, ax- 
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FIG. 10: Simulations with an externally-applied, axisymmet- 
ric magnetic field, (a) Kinetic and magnetic energy densities 
for a Re=116 (laminar), Rm= 70 (subcritical) simulation, (b) 
the resulting flow, (c) Kinematic and magnetic energy densi- 
ties for an Re=1803 (turbulent), Rm= 108 (subcritical) simu- 
lation, (d) The axisymmetric, time averaged velocity field for 
the turbulent simulation. The time average is over the time 
interval 0.3-2.4 t^, which is roughly 30 decorrelation times. 



isymmetric quantities and fluctuating quantities, 
B = (B) + and V = (V) + v, 



(30) 



where the brackets denote a time average over several re- 
sistive times. In practice, (B) and (V) are axisymmetric 
for sufficiently long time averages. Using these defini- 
tions, the time-averaged magnetic fields can be computed 
for laminar and turbulent flows, shown in Fig. 1111 

Both laminar and turbulent flows demonstrate toroidal 
fleld production and expulsion of poloidal flux. Laminar 
and turbulent results differ in several important ways, 
however, which are attributable to the currents being 
driven by MHD fluctuations. First, the toroidal field 
is greatly reduced in the turbulent run. The induced 
toroidal fleld is 6 times the applied field strength in the 
laminar flow and is only twice the applied field in the 
turbulent case. Second, the peak poloidal fleld is halved 
in the turbulent run, as shown in Fig. lllf a). Third, there 
is a net magnetic dipole moment associated with the 
induced fleld which is not present in the laminar case. 
These differences are partially the result of a difference 
between the mean flows in the two cases, but are mostly a 
due to a strong influence of the turbulence on the current 
generation. This can be interpreted in the context of a 
modiflcation to the mean-field Ohm's law, i.e. turbulence- 
generated currents are modifying the large-scale, mean 
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FIG. 11: (a) The magnetic field for a laminar flow described 
in Fig. 1101 The resulting total magnetic field (the sum of the 
externally applied field and those generated by the currents 
in the liquid metal) is shown as a multiple of the applied mag- 
netic field, (b) The time-averaged magnetic fields, scaled to 
the applied field, for a turbulent flow [Rm=107 (subcritical), 
Re=1803]. The peak internal poloidal magnetic fleld is 9.3 
times larger than the applied field. 
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FIG. 12: The wavenumber spectrum computed from fre- 
quency spectrum of fluctuations from Qt^ of flow (fltted with 
the red k~^^^ curve) output at a position (r ~ 0.75 a, 6 ^-^ 7r/2, 
4>=0) with a weak applied magnetic fleld of Bo '~ 51.3 G, fit- 
ted with the blue k~^^''' at low k and fc"^^/^ at large k. For 
this simulation, Rm — 130, and Re = 1450 and fluctuations 
are assumed to be due to convection of spatial variations in 
the field. The dispersion relation is a; = fc (V) . 



magnetic field. 

A turbulent EMF is possible because of the flow fluctu- 
ations and the magnetic fleld generated by the passive ad- 
vection of the applied magnetic fleld by the Kolmogorov- 
like turbulence in the velocity field. Fig. [12] shows the 
wavenumber spectrum as estimated from the frequency 
spectrum of the fluctuations in both V and i? at a fixed 
point in the simulation using the Taylor hypothesis to 
map frequency fluctuations to wavenumber to ~ k{V). 
It is clear that both the velocity field and magnetic fleld 
have an inertial range {k~^^^) and a dissipation scale, 
although the dissipation scales are at different values of 
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FIG. 13: The magnetic field (a) generated by the mean-fiow 
EMF (V) X (B). The peak internal poloidal magnetic field is 
11.2 times larger than the applied field. The magnetic field (b) 

generated by the turbulent EMF x (c) The currents 

associated with (a), and (d) the currents associated with (b). 
The turbulent flow has [Rm = 107 (subcritical), Re = 1803] 
and an externally generated magnetic field applied along the 
symmetry axis. The time average is taken over 2.lTa. Mag- 
netic fields are scaled to the strength of the applied field. 
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FIG. 14: The turbulent and helical fluctuations for the sim- 
ulation described in Fig. [T^ (a) The average of the squared 
turbulent fluctuations as a multiple of the peak squared mean 
flow, (b) The time average of the kinetic helicity fluctuations 
(v ■ V X v) as a multiple of the volume-averaged helicity of 
the mean flow. 



k. The k~^/^ scaling of the velocity field (the inertial 
range) is expected from the Kolmogorov theory of hydro- 
dynamic turbulence. The dissipation scale for the fluid 
turbulence is expected to be at ky ^ Re~^^'^ = 235 which 
is roughly the position of the viscous cutoff shown in Fig. 
1121 The limited inertial range at low k, is primarily due 
to constraints on long-time averages of the data imposed 
by computational speed. 

The k~^^^ scaling of the magnetic field corresponds to 



the weak-field approximation in which the induced mag- 
netic fluctuations are due to advection of the mean mag- 
netic field by the mean flow [s^ for k < k„ ^ Rni/a. 
The k~^^/^ power law results from a balance between 
the mean magnetic field advected by turbulence and the 
resistive dissipation of magnetic fluctuations. The dissi- 
pation scales are evident from the knee in the wave num- 
ber spectra of Fig.[T21 The spectrum is constructed from 
the power spectrum of the value of Br near the equa- 
tor. Consequently, the magnetic field gains structure at 
smaller scales as Rm increases, down to scale sizes of 
£cr = 2ti /ka- = 4.8 cm at Rm = 130. 

The simultaneous fluctuating magnetic and velocity 
fields can potentially drive current in a mean-field sense. 
The motional EMF can be written as 

V X B = (V) X (B) + (V) X b + V X (B) + V X b, (31) 

where the mean-fields have been separated from the fluc- 
tuating parts. The time averages must be taken over 
times long compared to a turbulent decorrelation time 
and comparable to the resistive diffusion time. Since the 
turbulent decorrelation time, tc ^ 0.05t^, integrating 
the induction term over several resistive times yields 

(v X B) = (V) X (B) + X b^ (32) 

An important question is whether the currents generated 
in the simulation are primarily due to the motional EMF 
associated with the mean-flow and the mean magnetic 
field, (V) X (B), or if there are also currents driven by 

the turbulent EMF <^v x b^ This can be investigated 

by examining the various terms in Ohm's Law 

E = r;J - (V) X (B) + (vxh^ (33) 

It is clear that in steady-state there can be no inductive 
electric field in the toroidal direction since the poloidal 
flux is constant. Axisymmetry precludes an electrostatic 
potential from driving current in the toroidal direction, 
and so the toroidal current can only be generated by the 
mean-flow and the turbulent EMF. Thus any currents 
driven in the toroidal direction contribute to the poloidal 
magnetic field. Fig. [TSFc) shows the currents driven by 
these fluctuations and their corresponding magnetic field 
(a). The fluctuation-induced magnetic field is 3.5 times 
larger than the applied field and comprises a third of the 
total field strength. 

It has been recently shown that an axisymmetric flow 
and axial magnetic field cannot induce a dipole moment 
in any simply-connected bounded system [43|. This is 
essentially due to the fact that the flow outside the con- 
ducting region is zero, while the streamlines of flow per- 
pendicular to the magnetic flux are closed and bounded 
within the conducting region. Only a turbulent EMF can 
create the dipole moment. With a weak applied field in 
a turbulent fluid, averaging over several eddy turnover 
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times and averaging along <f) eliminates the nonaxisym- 
metric component of the current, therefore the only non- 
trivial component of the dipole moment is 

^iz = <j) (fxrsm9J^. (34) 

The toroidal current generated by (^v x from Eq. [32l 

is shown in Fig. 113r d) and the associated dipole moment 
(antiparallel to the external field) is clearly seen in Fig. 
ITSTb). Alternatively, the EMF due to V and B gives rise 
to a hexapole magnetic field (in Fig. fTST a)). The resulting 
poloidal field reduces the surface magnetic field by 20%. 
The largest values of the turbulent toroidal current occur 
where the omega effect is also large. 

The EMF which generates the toroidal current asso- 
ciated with the induced dipole moment may very well 
resemble the currents driven by the well-known a-effect. 
The omega effect generates a toroidal magnetic field 
which would in turn would support a current of the form 
Jtf, = aB^. It is impossible to uniquely identify the cur- 
rent this way, however. A non-uniform /3-effect (change 
in local resistivity) could equally well explain the results. 
To do this would require separating the currents asso- 
ciated with the helical fluctuations from the non-helical 
fluctuations and this has not yet been done. A local anal- 
ysis of the turbulent helicity content in Fig. [HTb) shows 
that helical fluctuations exist that might be expected to 
drive a current through the a-effect. 

To study Ohm's law in the poloidal direction requires a 
full treatment of the poloidal electric field since an elec- 
trostatic potential is not ruled out by symmetry argu- 
ments. In MHD, the electrostatic potential is assumed 
to instantaneously adjust itself to ensure that V • J = 0. 
This can only be assured if the divergence of the motional 
EMFs is balanced by a spatially varying electric field 

V • E = -V^-S - = V • ((V) X (B) + X b^), 

(35) 

where $ and $ are electrostatic potentials due to the 
stationary EMF, and turbulent EMF respectively. Thus, 
a poloidal current can be associated with the mean-flows 
and the turbulent EMF respectively: 

Jpo; = ^ (-V* + (V) X (B)p,,) (36) 
Jpoi = <y (-V$ + (v X b^) . (37) 

When analyzing Ohm's law in the poloidal direction, it 
is necessary to first compute these potentials, which has 
been done for the poloidal currents in Fig. [131 

The simulations indicate there is a strong poloidal cur- 
rent, shown in Fig. [T3f d). associated with the fluctua- 
tions. The current acts to greatly reduce the toroidal 
magnetic field generated by a comparable laminar fiow, 
thereby reducing the toroidal field in the core. This re- 
sembles the diamagnetic 7-effect \A^, due to gradients 



in the turbulence intensity. The a-effect is the diago- 
nal part of a mean-field tensor: J — a a ■ B. The off- 
diagonal terms can also be written so that J = 0-7 x B. 
Figure [Hf a) shows the squared velocity fluctuations de- 
crease away from the axis of symmetry with the polar 
radius, p. For isotropic turbulence, the inhomogeneity in 
the fluctuations would give rise to a 7-effect of the form 
—(j-f{p) X Bt with 7 cx Vv^. The poloidal current due 
to turbulent diamagnetism would counteract the toroidal 
magnetic field caused by the omega effect. Comparison 
between Fig fT^f a) with FigfHTb) shows that regions of 
steep gradients in the turbulent fiuctuations correspond 
to regions of strong fluctuation induced poloidal current. 



VI. SUMMARY 

The role of turbulence in generating current and mod- 
erating the growth of magnetic flelds has been studied 
for the Madison Dynamo experiment using 3D numer- 
ical simulations. A simple forcing term has been used 
to model impellors in the experiment; at sufficient forc- 
ing the flow becomes turbulent. Two regimes were ex- 
plored: one with an external applied magnetic fleld and 
flow subcritical to the dynamo instability and one with 
no external fleld and super-critical flow. The role of the 
turbulence on current generation and self-excitation is 
marked. 

The onset conditions for the dynamo instability are 
governed not only by Rm but also by the magnetic 
Prandtl number Pm. At Pm ~ 1, the transition and sat- 
uration agree with laminar predictions and are considered 
laminar dynamos. At lower Pm, Rnicrit increases, con- 
sistent with a reduction in conductivity due to turbulent 
fluctuations. However, at higher Rm the character of the 
dynamo changes; its symmetry suggests that turbulence 
driven currents are important in the self-excitation pro- 
cess. The Pm values in the simulations are still orders of 
magnitude larger than in liquid-metal experiments (and 
for geo and solar dynamos) due to memory and speed 
limitations of computers, and so experiment support is 
critical for verifying these results. 

To quantify currents driven by fluctuations in the ex- 
periment, simulations of subcritical flows have been per- 
formed, and the currents driven by the turbulent fluctu- 
ations have been observed directly. The main effect of 
the turbulence on an externally-applied magnetic fleld is 
the reduction of fleld strength compared with those com- 
puted for laminar flows. The laminar two-vortex flow 
compresses the applied poloidal magnetic flux near the 
axis of symmetry and through toroidal flow shear cre- 
ates a strong toroidal magnetic fleld. Both effects are 
reduced in turbulent flows. The mean flow produced at 
large Reynolds numbers differs from its laminar counter- 
part, which accounts for some of the discrepancy between 
the build-up of toroidal fleld and flux compression of the 
poloidal fleld observed in the laminar and turbulent flu- 
ids. However, it has also been shown that a fluctuation 
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driven EMF drives current which modifies the large-scale We would like to thank C. Sovinec, S. Prager, J. 
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